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Abstract: Examples of joint probability distributions are studied in terms of Tsallis' nonextensive statistics both for correlated and 
uncorrelated variables, in particular it is explicitely shown how correlations in the system can make Tsallis entropy additive and 
that the effective nonextensivity parameter fjjy decreases towards unity when the number of variables N increases. We demonstrate 
that Tsallis distribution of energies of particles in a system leads in natural way to the Negative Binomial multiplicity distribution 
in this system. 
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1 Introduction 



It is very well known fact that whenever in a single variable exponential distri- 
bution 

/(,t) = i . exp (-^) (1) 



parameter A fluctuates according to gamma distribution, i.e., A ^ A' = A'(e) = 
^, where e is distributed according to 
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one obtains as result the following power-like distribution [T1I2| : 
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called Tsallis distribution and characterized by parameter q (q G (1,2), 
for g ^ 1 eq. (01) becomes the usual exponential distribution given by eq. 
dH)). Actually fluctuations can be also described by more general distributions 
(which induce wide spectrum of the so called superstatistics [T]), in all cases 
parameter q reflects the amount of fluctuations and is connected with their 
measure given by cj = Var{£) /(e)^. In the case of eqs. (EBSI) one has (e) = 2 — q, 
Var{£) = {q — and 
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This is result for the so called type B superstatistics |T]. In type A superstatis- 
tics, not accounting for A-dependent normalization f2], one gets q = 1 + u. In 
this case large fluctuations were corresponding to large whereas eq. (jlj) de- 
scribes all fluctuations using only limited range of parameter ranging from 
a; = for g = 1, up to cj ^ oo for g reaching its maximal allowed value of 
q ^ 2 (notice that for small values of to both approach lead practically to the 
same result: g ~ 1 + u). 

Such distributions are widely used to characterize systems with stochastic pro- 
cesses and are often associated with the existence of long-range correlations, 
with memory effects and with nontrivial (multi)fractal phase space structure 
PHj . Especially important factor in getting such distributions are all possible in- 
trinsic fluctuations which exist in the system under consideration P] ^ . In this 
paper we shall study examples of formal interrelations between fluctuations, 
correlations and nonextensivity, which can be of interest in physical applica- 
tions, especially in the field of high energy multiparticle production processes 
(cf., lEMSmi and references therein for other dynamical motivations). Notion of 
nonextensivity reflects the fact that Tsallis distribution (EI) has also its ori- 
gin in the so called nonextensive statistical mechanics (or information theory) 

* In fact in recent work (5j it was argued that Tsallis distribution does not imply dynamics with 
correlated signals but rather with signals occurring in nonstationary intervals of time. 
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based on Tsallis entropy, which depends on parameter q (becoming for q ^ 1 
the usual Boltzmann-Gibbs-Shannon entropy). It is normally nonextensive by 
amount proportional to q — 1, therefore q is named nonextensivity parameter 
P,4J. Anticipating some applications we shall in what follows have in mind dis- 
tributions of particles or their energies rather than some unspecified variables. 
Our goal is limited, the more general (but also very specialized) discussions on 
the role of correlations in obtaining Tsallis distributions can be found in jB], 
whereas in |[7||H1IS] one can find further discussions dealing with both fluctua- 
tions and correlations (in some specific scale-free approach) showing, among 
other things, that they can make Tsallis entropy additive. In next Section we 
show, using simple examples, that fluctuations in composite (but uncorrelated) 
systems can induce correlations and g 7^ 1, i.e., nonextensivity. On the other 
hand fluctuations in correlated system can make it apparently uncorrelated and 
extensive (see also Appendix A). When applied to multiplicity distributions of 
particles produced in collision processes of all kinds it is shown that they con- 
vert the usual Poissonian distribution to the co called Negative Binomial ones 
found in all high energy processes fTDj - this is shown in Section 3. The imme- 
diate practical applications of our results to some recent experimental data are 
presented in Appendix B. 

2 Tsallis distributions: fluctuations vs correlations 

2.1 Two random variable case 

To introduce and discuss correlations one has to deal with at least two particles. 
Let X and y denote therefore two independent random variables, each following 
its own exponential distribution given by eq. (P) and let their joint probability 
distribution be given by /(x, y) = f{x) ■ f{y). The corresponding Tsallis distri- 
bution can be obtained either by fluctuating the parameter A for each variable 
separately (in which case one obtains joint Tsallis distribution for uncorrelated 
random variables) or by fluctuating parameter A jointly for both variables (in 
which case one gets Tsallis distribution for correlated random variables) [IJ. It 
should be stressed that fluctuations lead always to Tsallis distribution, irre- 
spectively of the presence or absence of correlations (their introduction does 
not change the statistics). We shall now discuss in more detail the uncorrelated 
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and correlated cases separately. 



2.1.1 Uncorrelated random variables 



In this case we fluctuate independently parameter A in single particle probabil- 
ity distributions f{x) and f{y) as given by eq. dH), obtain in this way Tsallis 
distributions h{x) and h{y), and finally the joint Tsallis probability distribution 
in the form of: 



h{x, y) = h{x) ■ h{y) = 
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(5) 



However, such distribution of uncorrelated variables can be also obtained start- 
ing from joint distribution of two correlated variables x and y and introducing 
to it suitable fiuctuations in the way described before. Let us take, for example, 
the following two variable distribution. 



f{x,y) = 
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where {Q — 1) > describes (negative) correlations between variables x and y 
present for Q 7^ 1. It is characterized by the correlation coefficient p: 



Cov{x,y) ^ {xy) - {x){y) 
^JV ar{x)\JV ar{y) ^JV ar{x) ^JV ar{y) 
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(we have used here the exponential integral function Ei{z) = Ji° ^-^dt and 
incomplete gamma function r(0, z) = J^^dt). Fig. 1 shows how p depends 
on Q. Let us now fiuctuate in eq. (EI) parameter A' in the same way as in eq. 
(EI) with q being measure of these fiuctuations. It is straightforward to notice 
that one introduces in this way positive correlations between variables x and y 
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Fig. 1. Illustration of dependence of correlation parameter p as obtained by Eq, (jSJ on the nonexten- 
sivity parameter Q. 

and that for some strength of these positive correlations given by the condition 
q = Q they cancel the negative correlations introduced in eq. (EI) and we get 
uncorrelated variable distribution h{x,y) as given by eq. (0). This example 
shows that effects of correlations and fluctuations can cancel each other in final 
result . Let us notice at this point that correlations assumed in OBI) which lead 
to (ED were recently widely used (cf. jTT] and references therein) as example of 
the nonextensive rule for addition of energies resulting in Tsallis distribution. 
We argue that this is not true, namely that energies add always additively 
and formula used in [llj {n{x^y) = x + y + axy) leads to the formula for 
joint probability distribution, as our eqs. (ED or (EI), which does not provide 
distribution of the sum of energies (cf. ||T2] and references therein for different 
mappings leading to different composition rules). 

2.1.2 Correlated random variables 

To get correlated random variables one fluctuates parameter A in joint prob- 
ability distribution of two variables, {x, ?/}, as given by f{x,y) = f{x) ■ f{y). 

**In Appendix A we provide another example of this kind when similar approach results also in 
extensivity of the final g-entropies of the original distributions. It should be stressed that here q > 1 
whereas in the example presented in Appendix A we have always q < 1. 



5 



Performing it in the same way as in OHj) one gets the following joint Tsallis 
probability distribution: 



hix,y) = Jd£\^jj exp 
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with 
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(9) 



Notice that marginal probability distributions h{x) = / y)dy and h{y) = 
/ h{x, y)dx have in this case also form of Tsallis distributions but with noticeably 
difference in the exponent, which lacks now the factor q present in (jOI), i.e., they 
are identical with eq. (Ej) where exponent is by unity higher than in (jHI) (being 
equal 1/(1 — q) = q/{^ — q) + 1). The corresponding mean value and variance 
of the marginal distributions are equal to 

\ \'^(2 — n) 

(x) = (y) = and Var(x) = Var(y) = -ttt r, (10) 

^ ^ 3-2q ^ ^ (3-2g)2(4-3g)' ^ ' 



whereas the corresponding covariance and correlation coefficient corresponding 
to this joint probability distribution are equal to 

^""'"'^^^ (3 -^2^)^(4 -3,) ^""i = 1^ = 2^-1- (") 

It means then that correlation coefficient p is entirely given by the parameter 
which defines fluctuation of the random variable e (actually, in this case it is 
just equal to the relative variance of e given by cj, cf. eq. dH)). Notice that we 
can now express the nonextensivity parameter q via the correlation coefficient: 

^=1 + ^. (12) 

P+l 



In this way the correlation coefficient defines exponent in the Tsallis distribu- 
tion. However, in any possible applications one has to remember that eq. (1121) 
has been obtained only for correlations caused by fluctuations and therefore 
cannot be used for estimations of the role of correlations as such in Tsallis dis- 
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tributions. In Appendix A we present simple example of correlations making 
Tsallis entropy additive (cf. |71IH|IS|)- 



2.2 N random variables case 



It is straightforward to proceed to general case of N random variables {xi jsj} 
by fluctuating common parameter A in the corresponding initial exponential 
distribution /({xi^...^iv}) = f{xi) {with f{xi) being given by eq. dH)) following 
prescription given by eq. (j2I). We get in this case following joint probability 
distribution: 
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(13) 



where 



CN = ^m{i-2)q 

^ i=l 
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(14) 



It is straightforward to check that this distribution leads to marginal distribu- 
tions in form of eq. (01) for each of variables considered. 

Introducing A^-particle nonextensivity parameter one can formally rewrite 
eq. (1131) as: 
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= 1 - iV or qM=l+ , , ^ 1. (16) 
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Notice that, irrespectively of how large are single variable fluctuations (repre- 
sented by gi), they disappear in the multi-component systems with very large 
number of components TV * ^ ^ . 



3 Application to multiplicity distributions 



3.1 Boltzmann distribution and Poisson multiplicity distribution 



Suppose that in some physical process one has N independently produced sec- 
ondaries with energies {£"1 . . ^v}, each distributed according to Boltzmann dis- 
tribution, i.e., according to eq. dH) with x = Ei and A = {Ei). The corresponding 
joint probability distribution is then given by: 



/({^i,...,7v}) = -^-exp 



1 iV \ 

■^E^.i = n 



i=l 



1 



(17) 



For independent {£'i=i,...,iv} the sum E = E^i Ei is then distributed according 
to gamma distribution, 



gN{E) = 



1 



(18) 



with distribuant equal to 



N-l 1 



Ey-'^ , E 

■ expi 



A 



(19) 



* * * In jni) in thermodynamical context and where (7 < 1, one has — = |(1 — instead. 
Notice also that if (gi — 1) (N — 1) >> 1 one has approximately that 1/ (q^at — 1) = • 1/ (gi — 1), 
which seems to coincide with the notion of the extensivity of parameter ^ = l/{q — 1) discussed 
in |14ll3j . Notice also that (|15|) is multivariable distribution rather then distribution of the sum of 
variables discussed recently in the context of the conjectured g-central limit theorem [T^ . 
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We look now for such N that E^LqE^ < E < E,. Their distribution has 

known Poissonian form (notice that E/X= {N)): 

PiN) = Gj,+,{E) - Gn{E) = ^ ■ exp{-aE) = ^-J^ ■ exp(-(iV)).(20) 



In other words, whenever we have variables Ei, £"2, . . . , E^, ^'tv+i, • • • taken from 
the exponential distribution / [Ei) and whenever these variables satisfy the con- 
dition J^iLo Ei < E < EilV then the corresponding multiplicity has Pois- 
sonian distribution (actually this is precisely the method of generating Poisson 
distribution in the numerical Monte-Carlo codes). 



3.2 Tsallis distribution and Negative Binomial multiplicity distribution 



Suppose now that in another process one has again particles with energies 
{^i,...,iv} but this time distributed according to Tsallis distribution as given by 
eq.(IIBI) (therefore, according to our previous discussion they cannot be inde- 
pendent but are correlated in some specific way), 

(21) 
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h[{Ei iv})=C 
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with normalization constant Cn given by eq. (1141) ). It means that, according to 
our reasoning behind eq. ( |T3I) . there are some intrinsic (so far unspecified but 
summarily characterized by the parameter q) fluctuations present in the system 
under consideration. Because variables {£"^=1,...,^^} occur in the form of the sum, 
E = YliLi Ei, one can perform sequentially integrations of the joint probability 
distribution ( OTI) and, noting that 



hN{E) = %_i(E) 
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arrive at formula corresponding to the previous eq. (1181). namely 
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with distribuant given by 
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As before, for energies E satisfying condition J^^qEi < E < T.fJo' Ei, the 
corresponding multiphcity distribution is equal to 



P{N) = Hjv+i{E) - Hn{E) 



(25) 



and is given by the so called Negative Binomial distribution (NBD) (widely 
encountered in analyzes of high energy multiparticle production data of all 
kinds [IDl): 



P{N) = 
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where the mean multiplicity and variance are, respectively. 



{N) = 
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Var{N) 
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(26) 



{N) + {N)'.iq-l). (27) 



(For different way of deriving of NBD by using Tsalhs statistics see flB])- It is 
defined by the parameter k equal to: 

fc = ^. (28) 



Notice that for g ^ 1 one has k ^ oo and P{N) becomes Poisson distribution 
whereas for q ^ 2 one has k ^ 1 and we are obtaining geometrical distribution. 
As we have notice before, fluctuations described by parameter q result also in 
specific correlations described by parameter p given by eq. ( ITTl ). It means that 
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parameter k in NBD can be also expressed by the correlation coefficient p for 
the two-particle energy correlations (resulting from intrinsic fluctuations in the 
system), namely 

k = (29) 
P 



In should be stressed at this point that result ( OHl) coincides with our previous 
results in [T7] where we have already obtained NBD from fluctuations of the 
mean multiplicity in the Poisson distribution. It means that such fluctuations 
are equivalent to fluctuations leading to eq. ( |T^ which, following our reasoning 
presented in |2j , we would like to attribute to the fluctuations of temperature T 
for the whole system ^ * . We would like to close this Section with the remark 
that our result harmonizes with the known fact that whereas the generating 
function for Poisson distribution is exponential the corresponding one for the 
NBD has g-exponential form [THlfTB] (the relations between generating functions 
and probability distributions are the same as between P{Ei) and P{N) in our 
case, i.e., P{Ei) plays the role of generating function for distribution P{N)). 



4 Summary 



To summarize, let us stress that fluctuations that can be described by gamma 
distribution (or equivalent to it in the sense discussed in P where term of super- 
statistics has been coined for this purpose) lead always to Tsallis distribution. 
On the other hand, not every fluctuation results in correlation. This is true only 
for fluctuations of the whole multicomponent system {x + y) or (SjXi). Indepen- 
dent fluctuations of parameters in /(x, y) ~ [exp{—x/Xx)] ■ [exp{—y/Xy)] lead 
to distribution /i(x, y) = h{x)h{y) given by product of two Tsallis distributions 
with no correlations between variables x and y. Distributions obtained here 
differ from some many-particle distributions for composed systems {x^^i . . ^v} 
of the form exp(Eia:i) — ^ [1 + (1 — q) Y.iXi\^^l^^^'^^\ which occur as apparently 
natural (and simple) generalization of the observation that in single component 

•k-k-k^ To make this point more transparent let us notice that because (iV) = E/X therefore fluctuation 
of (A^) in Poisson distribution in is equivalent (for fixed E as in our case) to fluctuation of 1/A, 
i.e., in our case to fluctuation of (Ei). 
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systems fluctuations lead to the replacement exp(a:) — > [1 + (1 — q)x]^'^^'''^~'^^\ 
(and were also used to discuss correlations generated this way [TT]) but they 
do not lead to correct (marginal) single particle distributions. Finally, we have 
proved that energy correlations introduced in multiparticle system by fluctua- 
tions (which can be traced to fluctuations of temperature in this system, which 
is the place where energy is converted into observed particles in process known 
as hadronization) result in changing the corresponding multiplicity distribu- 
tions from Poisson to Negative Binomial ones and that this is equivalent to 
introducing fluctuations of the mean multiplicities in the Poisson distribution. 
The possible application of our flndings to some recent data on multiparticle 
production processes is presented in Appendix B. 

A Example of extensive Tsallis entropy 

In pj we have shown that nonextensivity leading to Tsallis statistics [i3j and 
characterized by parameter q can be caused by some intrinsic fluctuations ex- 
isting in the physical system under consideration. The corresponding Tsallis 
entropy is nonextensive. However, as discussed in jSlE IIHlEj (cf. also exten- 
sivity depends not only on the speciflc form of the entropy function used but also 
on the composition law according to which given composed system is formed 
out of its subsystem, i.e., on their possible correlations. In fact it is easy to 
demonstrate |Hr7IIH|PfT^ ) that by introducing to the physical system some spe- 
ciflc correlations (for example, correlations that are strictly or asymptotically 
scale invariant fTllH ipfT!]] ) one can make the corresponding entropy becoming ex- 
tensive. In what follows we shall illustrate this point by using as example simple 
gaussian probability distributions for single and two correlated variables. 

Let the single variable x distribution be of the gaussian form (with a = ^Var{x) 
being a parameter) 

Let us correlate this variable with another variable, using two variable gaus- 
sian distribution with correlations provided by parameter p as deflned in Eq. 
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(ED: 



f{x,y) 



exp 



— 2pxy + 
2(j2(l-/?2) 



(A.2) 



It is properly normalized, i.e., / / f{x,y)dxdy = 1, and has properly defined 
marginal probabilities, namely / f{x,y)dy = f{x). It is obvious that x and y 
are not independent because f{x,y) ^ f{x)f{y). 

The corresponding Shannon entropies for f{x) and f{x^y) are 



5', = ^[l + ln(27r) + 21n(a)], 



S^,y = 1 + ln(27r) + 2 ln{a) + In ^1 - /^^ 
respectively, being nonextensive by amount 



6S = S:,^y - 23:^ = ln^Jl - p^, 



(A.3) 
(A.4) 

(A.5) 



which depends on strength of correlation p (notice that 5S < 0) and becoming 
extensive only for uncorrelated system, i.e., when p = 0. 



The corresponding Tsallis entropies are equal to: 
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(A.6) 
(A.7) 

(A.8) 



Notice that for uncorrelated variables, i.e., for = 0, one gets in this case the 
usual result for Tsalhs entropy: 



5T={1- q)T 



(A.9) 
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However, it is obvious from eq. ( IA.8I) (cf. also Fig. lA.ll) that one can always 
find such value of correlation p for which 5T = and Tsallis entropy becomes 
extensive: 



9-1 



2(27r)^v^-(27r)''-ig 



1-9 



(A.IO) 



or, in crude approximation, 

— ^ q; < q < I. 



(A.ll) 



Notice that, according to eq. ( IA.5I) . both positive and negative correlations 
result in some (negative) imbalance of Shannon and Tsallis entropy. One has 




Fig. A.l. Left panel: \p\ versus q: full line - exact formula (|A.10|) : dotted line - approximate formula 
(|XTT1) . Right panel: y/l - versus q. Notice that < q < 1. 



to stress at this point that what we have shown here is only a kind of formal 
exercise presented for illustration use only. It does not prove that correlations 
of some type lead to Tsallis type distributions. Variables can be correlated in 
any way and this does not depend on the distribution. 



B Example of possible practical applications 



Results presented in this paper have practical application in the field of high 
energy multiparticle production reactions, especially those originated by col- 
lisions of heavy nuclei, which are of particular interest as potential source of 
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production of new state of matter, the so called Quark Gluon Plasma (QGP) 
(cf. references in POICT^^PB]). Recently number of works POICT^'SJ have 
demonstrated the existence in such reactions event-by- event fluctuations of the 
average transverse momenta {p) per event. The following quantities were consid- 
ered: Var{{p)) / {{p)Y Q'lid (ApiApj) / {{p))'^ . These quantities, as we advocate, 
are fully determined by u as defined by eq. (jll) ,which in our case translates into 
fluctuations of the temperature T of hadronizing system - a vital observable 
when searching for QGP ^ ^ * ^ ^ . The results obtained in can be 

interpreted as fluctuations of the temperature T of the hadronic matter being 
produced (these are fluctuations for the whole event or its part (cluster) but 
not for the particular (single) particles in an event). 

To show this let us consider the case of N^y events with Nk particles in the k^^ 
event. Introducing the following notation: 



Ck=J:J:iP^-m)■{PJ-m). m) 

* 3 

I N^v I Nk 

{{p)) = = -vr- E {p)k] where {p)k = ^ J2 Pt^ (B.2) 



Nev k Nk 



we have that 



Adding and subtracting the same term (p)l it can be written as 



1 Ne^ 

c= — E 



'eu k 



1 Nk Nk 

Y^T.J:p.n + {(p)l-(p)l)-iip)y 



(B.4) 



•k-k-k-k *Qgnerally speaking, analysis of transverse momenta px alone indicates very small fluctuations 
of T. On the other hand, as reported in the measured fluctuations of multiplicities of produced 
secondaries are large (i.e., multiplicity distributions are substantially broader than Poissonian). Our 
analysis of NBD applied to observed multiplicity distributions show that this can result in large 
fluctuations of T, cf. [25) . 
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If particles in the event are independent then 




(B.5) 



and we have that 




{{p)f)=Vari{p)), 



(B.6) 



or that 



c 



Var{{p)) 



Var{T) 



= UJ. 



(B.7) 



The above formulas can be checked against data obtained at Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven Nat. Lab where Au nuclei are impinging 
at each other with center of mass energy 200 GeV per nucleon. Data taken by 
STAR experiment |22j for centrality 30 - 40% {Npart ~ 100) give cj = 4 • lO'l 
The respective values for data obtained by other experiment, PHENIX, [23] 
are: 2.2 • 10~^, 2.4 • 10~^, 3.6 • 10~^ and 4.9 • lO"'^ for the respective centralities: 

- 5%, - 10%, 10 - 20% and 20 - 30%. The more detailed analysis of the 
RHIC data along the line presented here is, however, out of the scope of the 
present paper and will be presented elsewhere. 

Partial support of the Polish Ministry of Science and Higher Education (grants 
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